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Polynomially-large ground-state energy gaps are rare in many-body quantum systems, but useful 
for adiabatic quantum computing. We show analytically that the gap is generically polynomially- 
large for quadratic fermionic Hamiltonians. We then prove that adiabatic quantum computing can 

CI 

realize the ground states of Hamiltonians with certain random interactions, as well as the ground 
(3JT)' states of one, two, and three-dimensional fermionic interaction lattices, in polynomial time. Finally, 

r-{ , we use the Jordan- Wigner transformation and a related transformation for spin-3/2 particles to 

show that our results can be restated using spin operators in a surprisingly simple manner. A 
direct consequence is that the one-dimensional cluster state can be found in polynomial time using 
CC ' adiabatic quantum computing. 
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I. INTRODUCTION 



■ Adiabatic quantum computing (AQC) [1| is an approach to quantum computation where a problem is encoded 
as the ground state of some Hamiltonian Tip. It is assumed that it is feasible to prepare a physical system in the 
ground state of some simple Hamiltonian 7io, and then evolve the Hamiltonian slowly from Tio to Tip. Under the 
right conditions and if the evolution is done sufficiently slowly Q , then at the end of the evolution the state of the 
system is the ground state of Hp. Measurement of this final state reveals the solution to the original problem. As 

■ an approach to quantum computing, AQC is known to be equivalent to standard gated quantum computing, in that 
each can be efficiently simulated by the other [!, Hj]. Also, a simple Hamiltonian evolution corresponds to the Grover 

^^O ■ search algorithm [f| . 

AQC succeeds in polynomial time if the inverse of the ground-state energy gap is bounded by a polynomial in the 
. . problem size. A typical Hamiltonian must fit exponentially-many energy levels into a polynomial-sized energy range, 
so most energy gaps must be exponentially small. It is not clear a priori why the ground-state energy gap should 
ever be larger than the rest. Let us be more precise. Mathematically, we can construct a Hamiltonian with any given 
set of 2 n energy levels. Theorem lI.il (below) establishes that large ground-state energy gaps are rare among choices 
of energy levels. 



Theorem 1.1 (Large ground-state energy gaps are rare). Consider uniform random choices of energy levels for 
a Hermitian operator on n qubits, under the restriction that the ground-state energy is zero and the energies are 
^ , contained in the unit interval. The fraction of these choices with a ground-state energy gap greater than 2~™/ 2 tends 
■ ■ ' to e -2 for large n. 

Proof. The ground-state energy gap is larger than some e > provided the 2™ — 1 non-zero energy levels are selected 
from the interval (e, 1]. Thus the fraction of choices of energy levels with a ground-state energy gap of at least e is 
(1 — e) 2 " -1 . Let us choose e = 2~"/ 2 , then we have 



f 1 - 2-"/ 2 J = M - 2-"/ 2 J (1) 
-(e- 1 )^ (2) 

~ -2" /2 fo\ 
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□ 

In fact, since the dimension of the problem is exponentially large in the number of qubits, it is difficult to even 
determine the minimum ground-state energy gap for large problems. Therefore, it would be useful to identify a class 
of Hamiltonians meeting the following requirements: 

1. The class should have large ground-state energy gaps. 

2. The class should have many degrees of freedom and allow many interactions between qubits, so that it represents 
diverse problem- instances for AQC. 

3. The class should include a simple 7io for AQC, and be convex, so that Tt(s), where 

H(s) - (1 - s)H + sH P , (4) 

stays within the class for s £ [0, 1]. 

We show that quadratic fermionic interactions can be used to define a class of Hamiltonians meeting all three 
requirements. In Section[TTl we develop background on the fermionic commutation relations, necessary for Section [TTTT 
where we identify a class meeting all three requirements. We show that, provided the interaction coefficients are 
bounded, the fraction of Hamiltonians with a ground-state energy gap greater than 2x/n for any x > tends to e~ x 
for large n. For random Hamiltonians whose interaction coefficients are constructed from standard normal Gaussian 
distributions, we show the ground-state energy gap is 0(l/y/n). In Section HVl we analytically bound the ground-state 
energy gap for two and three dimensional lattices of interacting fermions. We then show that AQC can find the ground 
states of certain random Hamiltonians and these two and three dimensional fermionic lattices in polynomial time. 
Finally, in Section [V] we derive alternate representations of these Hamiltonians using the Jordan- Wigner transform 
and a related transform for spin-3/2 particles, and show that the one-dimensional cluster state may be obtained in 
polynomial time using adiabatic quantum computation. 



II. BACKGROUND ON THE FERMIONIC COMMUTATION RELATIONS 

For a detailed exposition on properties of the fermionic commutation relations (FCRs), see, e.g., @, @, Hi- Here we 
only highlight some essential points, mostly without proof, that is needed to develop the results in following sections. 
The FCRs on a set of linear operators {cj : j = 1 . . . n} are 

{ c 3>4} = S 3,k, {Cj,C k } = 0, (5) 

where the bracket notation indicates the anti-commutator {x, y} — xy + yx, and Sjk equals the identity operator if 
j = k and zero otherwise. The superscript dagger denotes the Hermitian adjoint. A consequence of the FCRs is that 
{cj, : j — 1 . . . n} are creation and annihilation operators that anticommute. 
Suppose we have a Hamiltonian of the form 

n 

H = Y,C j c]c j , (6) 
i=i 

where the coefficients Cj are positive and real. All the terms in H commute, and the j th term has eigenvalues and Cj. 
Now, take the sum of elements in each possible subset (including the empty set) of {Cj : j = 1 . . . n}. The 2 n resulting 
values are the eigenvalues of H. In particular, the ground-state energy of Ji is zero, and the ground-state energy gap is 
the least non-zero coefficient Cj . To decide whether an arbitrary value is an eigenvalue of Ti for arbitrary coefficients 
is NP-complete however, as it is equivalent to the subset-sum problem (also known as the knapsack problem) 

We can write many Hamiltonians in the form of © using Theorem III. II (below) , originally due to Lieb et al. 0] . 
Suppose we have a quadratic fermionic Hamiltonian Ti, defined as 

n 

n = A ^ { c l ck - c i c i) + B ^ k { c H - c i ck ) > ( ? ) 

j,k=l 

for some set of real coefficients Aj.k and Bj^- For convenience, we gather the coefficients Aj^ and Bj^ into real 
n x n matrices that we label A and B. If B = 0, then Ti represents a Hubbard model with no on-site interactions, for 
instance electrons in metals [Io| or graphene 11 1. If A and B are tridiagonal, then Ti represents a one-dimensional 
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chain of interacting spin-1/2 particles. If A and B have three non-zero super- and sub-diagonals, then TL represents 
a chain of interacting spin-3/2 particles [l2||. Also, we can see using the FCRs that different choices of A and B may 
represent the same Hamiltonian. In particular, for any given Hamiltonian, A can be chosen to be symmetric and B 
anti-symmetric. 

Theorem III . 1 1 establishes that we can write (O in the form of ([6]), added to a multiple of the identity. Thus we can 
easily find the first few eigenvalues of TL, and in particular its ground-state energy gap. 

Theorem II. 1 (Lieb et al., 1961). Consider a quadratic fermionic Hamiltonian as in where A is an n x n 
real symmetric matrix, B is an n x n real anti- symmetric matrix, and the operators {ck : k = 1, . . . , n} satisfy the 
FCRs. Then we can find A 2 diagonal and X unitary so that X(A — B)(A + B) = A 2 X , and Y unitary so that 
Y(A + B)(A - B) = A 2 Y. Define the operators {ly : j = 1, . . . ,n} by 



i = 1 ELi ( X 3k + Y ik ) c k + (X jk Y jk ) c\ 
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l) = I ELi {Xjk ~ Y jk ) c k + (X jk + Y jk ) c{ 



J 1 • (8) 



Then {r)j : j = 1, . . . ,n} satisfy the FCRs, and 



H = J2 2 *ivfa- (E^j^: (9) 

where Aj denotes the j th entry on the diagonal of the matrix A and /2 71 is the identity operator. 

Proof. See the appendix. □ 

Theorem III. II was used initially by Lieb et al. to find the spectrum of the one-dimensional XY model, and sub- 
sequently has been used, for instance, in the analysis of the one-dimensional model of free electron transport 0]. 
Quadratic fermionic Hamiltonians as in ([7]) have also sparked recent interest because of their application to quantum 
complexity theory. If one takes a set of gates defined by U — cxp(iTLt) for some t and a constant quadratic fermionic 
Hamiltonian Tt in the form of ([7]), then one obtains a set of gates that resembles a universal set, but in fact may be 
classically simulated [lH . Broader sets of gates that can be classically simulated have been identified [3, Ufl . To 
classically simulate an evolving Hamiltonian, it has been shown (3) that Hamiltonian evolutions may be efficiently 
approximated by discretizing the evolution into a sequence of short, constant Hamiltonians. Theorem III. II has also 
been used to find efficient sets of quantum gates for computing properties of quadratic fermionic Hamiltonians in the 
form ([7]) [la ]. Further, the relationship between a vanishing energy gap and discontinuity in the ground state has 
been studied for these Hamiltonians fl71 |. 

The Jordan- Wigner transformation applied to the Hamiltonian evolution 

n n — 1 

H(s) = (1 - s) £ a] + a £ o*o* +1 (10) 
j'=i i=i 

transforms 7i(s) into a quadratic fermionic Hamiltonian in the form of ((7]), thus providing a means for determining 
the spectrum for any s € [0,1]. In fact, this evolution exhibits a second-order quantum phase transition, and a 
ground-state energy gap that decreases as 0(l/n) [TH, [l||, which is surprisingly large. Yet much broader classes of 
quadratic fermionic Hamiltonians in the form of ([JJ also have remarkably large ground-state energy gaps. 

Since (A + B) l[ = A — B, and since the singular values of a matrix M are the square roots of the eigenvalues of 
M^M, we see that Aj from is a singular value of A + B. Further, Cj in ^ can be defined to be 2Aj. Thus, if A + B 
is non-singular, then twice the least singular value is the ground-state energy gap of Ti.. If A + B is singular, then 
Tt has a degenerate ground state, and the least non-zero singular value is the energy gap between the ground-state 
subspace and the higher energy levels of the Hamiltonian. In any case, since A + B has only n eigenvalues, in contrast 
to TL which has 2™, we might expect that often the least singular value of A + B is not exponentially small in n. Then 
the ground-state energy gap of TL would not be exponentially small. In the next section we state and prove more 
precise formulations of this claim. 



III. A CLASS OF HAMILTONIANS MEETING THE REQUIREMENTS OF AQC 



The class of quadratic fermionic Hamiltonians represented by ([7]) in fact meets all three requirements to be useful for 
adiabatic quantum computing. To establish that the ground-state energy gaps are large (Requirement [1]), we provide 
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two theorems. In Theorem IIII.ll (below), we take a particular distribution of coefficient matrices A and B under 
the restriction \\A + B\\ 2 < 1, and establish that the ground-state energy gap is 0(l/n). Here all matrix norms are 



assumed to be the two- norm, which is equal to the square root of the largest singular value. Then, in Theorem IIII. 31 
we show that the ground-state energy gap is 0(1/ 'y/n) for Gaussian-distributed interaction coefficients. 

Theorem III.l (Ground-state energy gaps of quadratic fermionic Hamiltonians with bounded coefficients). Choose 
a real diagonal n x n matrix £ uniformly at random with entries in the unit interval, and choose U and V according 
to any probability distribution over orthogonal n x n matrices. Then C = UYjV^ represents a distribution over all 
real matrices with \\C\\ 2 < 1. Take A to be the symmetric part of C and B to be the anti- symmetric part of C, e.g. 
A = (C + C^)/2 and B = (C — C^)/2, and let Tt be defined as in The probability that the ground-state energy 
gap of Ti is greater than 2x/n, for any x > 0, tends to e~ x for large n. 

Proof. If the ground-state energy gap Ti. is greater than 2x/n, then the singular values of C are contained in the 
interval (x/n, 1). The fraction of choices for £ where this is true is 







)•- 









(11) 



which tends to e 1 for large n. □ 

To determine the ground-state energy gap for Gaussian-distributed interaction coefficients, we first need the fol- 
lowing theorem about random matrices, due to Edelman pol . Corollary 3.1]. 

Theorem III. 2 (Edelman, 1988). Let C be an n x n matrix, whose elements have independent Gaussian distributions 
with mean zero and unit variance. We denote such distributions as N(0,1). Let q be the least singular value of C. 
Then for large n, ns 2 converges in distribution to 

p(x) = ^e-(*/ 2 +^) . (12) 

Since n^ 2 has a probability distribution that is asymptotically independent of n, it follows that <r = 0(l/y/n). Also, 
(I12[) implies that s ^ with probability one. Similar results for other ensembles of random matrices are known [2lj . 
Let us now apply Theorem IIII. 21 to quadratic fermionic Hamiltonians in the form of (JT)) . 

Theorem III. 3 (Ground-state energy gaps of quadratic fermionic Hamiltonians with Gaussian coefficients). Let C be 
an nxn matrix with independent N( 0,1) coefficients, let A be the symmetric part of C , and let B be the anti- symmetric 
part of C , so 

(13) 

and C = A + B. Define 

n 

H = A ^ k { c \ Ck - c i c l) + B ^ ( C K ~ c j°k) > (I 4 ) 

j,k=l 

and let 7 be the ground-state energy gap ofH. Then, for large n, n~f 2 /4 converges in distribution to p(x) defined in 



Proof. By Theorem IIII. 21 if 7/2 is the least singular value of C, then nj 2 /4 converges in distribution to p(x) for large 
n. Theorem IIII . 2 1 also implies that C is non-singular with probability one, so 7 is the ground-state energy gap. □ 

Since rry 2 /4 has a probability distribution that is asymptotically independent of n, 7 = 0(l/y/n). Recalling 
Theorem 11.11 we see this is a remarkable property. Since there must be 2™ distinct energy levels in an energy range 
of C(?i 2 ), most of the energy gaps must be exponentially small. In fact it can be shown that the Hamiltonians in 
Theorem IIII.3I are with probability one non-degenerate, so these exponentially small gaps are also non-zero. Figure [1] 
illustrates the difference between the distribution of the ground-state energy gaps and the rest of the gaps for 1000 
randomly-generated 10-qubit Hamiltonians, and indeed the ground-state energy gaps are typically much larger than 
the other gaps. 

To establish that quadratic fermionic Hamiltonians as in (JT]) have many interaction degrees of freedom (Rcquirc- 
ment[5|), we observe that for n qubits there are n(n+ l)/2 continuous degrees of freedom in choosing a real symmetric 
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FIG. 1: The ground-state energy gap distribution is compared to the distribution for the other energy gaps. All the energy 
levels are computed for 1000 random n = 10 (10-qubit) Hamiltonians. Each Hamiltonian is chosen randomly as described in 
Theorem IIH. 31 As predicted by Theorem IIII. 31 the ground-state energy gaps are much larger than the other gaps. 



n x n matrix A, and n(n — l)/2 continuous degrees of freedom in choosing a real anti-symmetric matrix B. The total 
degrees of freedom are thus n 2 . 

Finally, it is obvious that if two quadratic fcrmionic Hamiltonians are in the form of (J7|), then their sum is a 
quadratic fermionic Hamiltonian in that form (Requirement [3]) . We choose, as our initial Hamiltonian for AQC, 

n n 

^o = E ( 2c h - = E { c h - c A) ■ ( 15 ) 

J'=l i=i 

The ground state of this Hamiltonian is easy to construct. For example, for electrons in a metal, Cj is the annihilation 
operator for electron occupation at site j, and then the ground state of Ho is the state with each site unoccupied. 
Then the whole Hamiltonian evolution H.(s) is a quadratic fermionic Hamiltonian in the form of 10, where 

n 

H(s) = (1 - s)H + s E [Aj,k (ctc fc - Cjcfy + Bj, k (c]c{ - c 3 c k )] . (16) 

j,k=l 

To find the ground-state energy gap of Tt(s), we define 

A(s) = (1 - S )J 2 n + sA , 

B(s) = sB . (17) 

Then we can rewrite Equation (fTT))) as 

n 

H(s) = J2 AM (ct Cfc - c 3 c\) + B hk {s) (c]c\ - c jCk ) , (18) 
3 -,fc=i 

and twice the least non-zero singular value of A{s) + B(s) is the ground-state energy gap of Tt(s). 

We cannot directly use Theorem IIII. 31 to establish that the ground-state energy gap is large for all of H(s) in (|16[) . 
since Theorem IIII.3I is a probabilistic result for a single random Hamiltonian. However, for several special classes of 
Hamiltonians in (|16|) , we can establish that the gap is large throughout an evolution. 



IV. HAMILTONIAN EVOLUTIONS WITH LARGE GROUND-STATE ENERGY GAPS 

There are several classes of A and B for which we can easily find the least singular value of A(s) + B(s) as defined 
in p7p . and thus find the minimum ground-state energy gap of the evolution in (|16p . We do this for certain random 
choices of A and B, and also two and three-dimensional interaction lattices. 



0.2 0.4 0.6 0.8 1 
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FIG. 2: Eigenvalues of "H(s) as a function of s, where "H(s) was defined as in Theorem HV. II with A = CC 1 /n for n = 8. 
The division by n is so that ||7i(0)[| ~ ||7i(l)||, resulting in a better visualization although making the ground-state energy 
gap C(l/n 2 ) instead of 0(l/n). We can see that the ground-state energy gap is linearly decreasing with s as predicted by 
Theorem lIV.ll and is much larger than most of the other energy gaps. 



If A + B is symmetric positive semi-definite, we can determine the least singular value of A(s) + B(s) from the least 
singular value of A + B. In that case B = and thus Tt(l) represents a Hubbard model with no on-site interactions. 
In order to choose random samples of symmetric positive semi-definite matrices, we first choose C to be an n x n 
matrix with independent random N(0,1) elements, as before, and then set A = CC^ and B = 0. 

Theorem IV. 1 (Ground-state energy gaps for random Hamiltonian evolutions). Let C be an n x n matrix with 
independent N( 0,1) coefficients, and let A = CC^ . Define 

n 

H(s) = (1 - s)Ho + s Y, A i,k (c]c k - Cj cl) . (19) 

j,k=l 

Then the ground-state energy gap for TL{s) is 

7 (s) - 2(1 - s) + s 7 , (20) 
where 717/2 converges in distribution for large n to p{x) defined in 171 



Proof. Let us label the least singular value of C as y/ r y/2. Then by Theorem IIII.21 for large n, 717/2 converges in 
distribution to the density function p(x), and is non-zero with probability one. Since A — CC^ is symmetric positive 
semi-definite, 7/2 is its least singular value. Thus 7 is the ground-state energy gap of H(l). 
Define 

A(s) = (1 - s)I n + sA , (21) 

then the ground-state energy gap j(s) of TC(s) is twice the least non-zero singular value of A(s). Notice A and / are 
symmetric positive semi-definite, and diagonal in the same basis. Since 7/2 is the least eigenvalue of A, we have 

^-<i-H.(i)- <-> 

□ 
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Figure [2] illustrates the energy levels for an 8-qubit instance of a random evolution as in (|19p . Evidently the 
ground-state energy gap is larger than the other gaps throughout the evolution. 

If A + B is not symmetric, then the eigenvalues of G(s) = (A(s) + B(s))(A(s) — B(s)) are of interest. We explore 
three classes of matrices where finding these eigenvalues is easy, namely, generalizations of the one-dimensional XY 
model, and generalizations of two and three-dimensional lattices of interacting fermions. 

The one-dimensional XY model is expressed with A and B matrices such as [H, p. 413]: 



A=l 

2 



V i 



\ 



1 

/ 



B = l 
2 



( ° 1 

-1 



1 

1 



V i 



\ 



-1 J 



(23) 



where omitted entries are zero. The ground state of the one-dimensional XY model can be found in polynomial time 
with AQC OS EH- 

In fact, this holds for more general choices of A and B. An essential property of the definitions in 
(|23p is that each row is a cyclic shift of the previous row. Such matrices are called circulant. For n qubits, it is easy 
to check that there are n degrees of freedom in choosing a symmetric circulant matrix A and anti-symmetric circulant 
matrix B. We show that AQC can find the ground state of any Hamiltonian with circulant A and B matrices in 
polynomial time. 



Theorem IV. 2 (Ground-state energy gaps for circulant A and B matrices). Let 



n 



E 

j,k=l 



44 



CjC k 



(24) 



where A is a real circulant n x n symmetric matrix, and B is a real circulant anti- symmetric n x n matrix. Then the 
ground-state energy gap of TL is bounded from below by a polynomial in 1/n. 

Proof. Circulant matrices form a commutative ring [12, p. 201], so if A and B are circulant, then so is G = (A + 
B)(A — B). Also, G is symmetric positive semi-definite, and the ground-state energy gap of 7i is twice the square 
root of the least non-zero eigenvalue of G. 

Circulant matrices also have the nice property that their eigenvalues are given by the discrete Fourier transform of 
their first column [23l . p. 124]. The n x n Fourier transform matrix is 



F„ 



( 



\ 



3 (0-0)27ri/n 
D (l-0)27ri/n 



e (0-l)27ri/n 
e (M)27rj/« 



3 (0-2)27ri/n 
,(X-2)2m/n 



e (2-0)27ri/« e (2-l)27ri/n e (2-2)27rz / n 



= ((n-l)l-0)27ri/n ((n-l)-l)27ri/n {(n-\)-2)2ixi j n 



\ 



/ 



(25) 



Recall we label the eigenvalues of G as A^. Labeling the entries in the first column of G as gf., we can write: 



/A?\ 




(9i\ 


A| 1 




92 


A| 


= F n 


93 






Kg'n) 



(26) 



Also, G is symmetric so gu 



9n+2-k for k > 2. Then we have, for n odd: 
1 



(n+l)/2 

3=2 
(n+l)/2 

3=2 
(n+l)/2 

9i + E 9 i cos 

J=2 



,((i-l).(fc-l))27Ti/n , p ((n-i+l)■(fe-l))27r^/n , 



-(0-l).(fc-l))27i-l/n 1 



27r(j-l)(fc-l) 



(27) 
(28) 
(29) 
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For n even we have a leftover term, but it simplifies: 

i 



.9i + .9 n /2+ie 



((n/2)-(fc-l))27ri/rc 



n/2 
3=2 



27T(i-l)(fc-l) 



n/2 

3l + (-l) (/S_1) 3n/2+l+5^ft 
J=2 



27r(j-l)(fc-l) 



(30) 



(31) 



For n > 1, we have l/n < 1/^/n. So whether n is even or odd, Taylor expansion of the cosine makes it clear that 
if Afc 7^ then is bounded from below by a polynomial in l/n. So the ground-state energy gap of H is bounded 
from below by a polynomial in l/n. □ 



We can extend this result to a whole Hamiltonian evolution. 
Corollary IV. 3 (Hamiltonian evolutions with circulant A and B matrices) 



Let 



H{s) = (1 - s)H 



- J2 [ A ' J - ( C J C * 



,t„t 
■j c k 



(32) 



where A is a real circulant n x n symmetric matrix, and B is a real circulant anti- symmetric n x n matrix. Then the 
ground-state energy gap of TL{s) is bounded from below by a polynomial in l/n and s. 



Proof. First rewrite H(s) as in (fT5|) . The elements of the matrix 

G(s) = (A( S ) + B(s)) 



B(s] 



(33) 



are quadratic functions of s, and A(s) and B(s) are circulant. Then we can follow the proof of Theorem IIV.2I to see 
that the eigenvalues of G(s) are bounded by a polynomial in s and l/n, and so the ground-state energy gap of TL{s) 
is bounded from below by a polynomial in l/n and s. □ 

Circulant coefficient matrices other than the XY model include scenarios such as non-nearest neighbor interactions 
on a one-dimensional chain of interacting fcrmions. The restriction that A and B are circulant imposes the requirement 
that the interaction strengths are independent of position, and it imposes periodic boundary conditions. It should be 
noted that while circulant matrices yield to elegant analysis, these results could be extended to interaction matrices 
derived from other boundary conditions. For instance, if G is symmetric Toeplitz tridiagonal (Tocplitz matrices are 
those with constant diagonals), then its eigenvalues may be found analytically [13, p. 158]. Analysis of the two and 
three-dimensional interaction grids then builds on the onc-dimensional analysis in exact analogy to the circulant case. 

Now let us consider the case of a two-dimensional lattice of interacting fermions with periodic boundary conditions. 
The A and B matrices then have a block structure such as 



A = 



f Aq I 
I A 



V i 



I 



I A I 
I A J 



B = 



/B I 
f -/ B Q 



\ I 



-I 



I 



-I B I 
-I B ) 



(34) 



where Bo and Aq are as in (|23|) . Evidently each block is circulant, and A and B arc circulant in the blocks. Such 
matrices are called "block circulant with circulant blocks" (BCCB). Let us assume that the blocks are pxp, and there 
are q blocks per row, so n = pq. 

Theorem IV. 4 (Ground-state energy gaps for BCCB A and B matrices). Let A and B be n x n BCCB matrices 
with pxp blocks and q blocks per row, so pq = n. Assume A is real symmetric and B is real anti- symmetric. Define 
the quadratic fermionic Hamiltonian TL with A and B as in ([7]). Then the ground-state energy gap of Tl is bounded 
from below by a polynomial in l/n. 

Proof. It is easy to check G = (A + B)(A — B) is BCCB, and since G is symmetric, it is block symmetric and its 
blocks are symmetric. Then T^GT is diagonal 25j, where 



T 



/ F p 




(F q \ 




y 


i "J 


V 







(35) 
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and P y is a permutation matrix that re-orders the columns as (1, p + 1, 2p+l, (g— l)p+l, 2, p + 2, 2p + 2, ...). 
The matrix T diagonalizes G by first diagonalizing the blocks of G, then reordering the rows and columns so that the 
matrix is block diagonal with circulant blocks, and finally diagonalizing those blocks. Geometrically, this procedure 
can be thought of as a Fourier transform first along the horizontal axis of the lattice and then along the vertical axis. 

The eigenvalues of G are then the eigenvalues of symmetric circulant matrices, whose entries are the eigenvalues of 
symmetric circulant matrices. Following the proof of Theorem IIV. 21 we can see that the non-zero eigenvalues of G 
are bounded from below by a polynomial in 1/n. □ 

We can extend this result to a whole Hamiltonian evolution: 

Corollary IV. 5 (Hamiltonian evolutions with BCCB A and B matrices). Let 

n 

H(s) = (1 - s)Ho + a J2 [ A i,k (''/-A " Cjcl) + Bj :k (c)c{ - Cj Ck)] , (36) 

where A is a real BCCB n x n symmetric matrix, and B is a real BCCB anti- symmetric n x n matrix. Then the 
ground-state energy gap of TL{s) is bounded from below by a polynomial in 1/n and s. 

Proof. The proof is analogous to the proof of Corollary IIV. 31 □ 



We can even do the three-dimensional lattice of interacting fermions. Then the A and B matrices are as in (|34j) . 
but A and B are BCCB instead of circulant. Let us call these matrices "block circulant with BCCB blocks", or 
(BC) 2 CB. 

Theorem IV. 6 (Ground-state energy gaps for (BC) 2 CB A and B matrices). Let A and B benxn (BC) 2 CB matrices. 
Let the number of blocks be r, and the BCCB blocks contain q circulant subblocks each p x p, so n = pqr. Define the 
quadratic fermionic Hamiltonian TL with A and B as in Qj. Then the ground-state energy gap ofH is bounded from 
below by a polynomial in 1/n. 

Proof. The proof is analogous to that of Theorem IIV. 41 but we set 

F P \p y \ F 1 \P Z \ F r 



where P z is a permutation matrix that re-orders the columns as (1, pq + 1, 2pq + 1, . . . , (r — l)pq + 1, 2, pq + 
2, 2pq + 2, . . . ). This transformation first diagonalizes the BCCB blocks, then permutes rows and columns to obtain 
a block-diagonal matrix with circulant blocks, and diagonalizes the remaining blocks. □ 

We can extend this result to a whole Hamiltonian evolution: 

Corollary IV. 7 (Hamiltonian evolutions with (BC) 2 CB A and B matrices). Let 

n 

H(s) = (1 - s)H + s J2 [Aj, k (ctefc - Cj c[.) + B jtk (c]c\ - c 3 c k )] , (38) 
j,k=l 

where A is a real (BC) 2 CB n x n symmetric matrix, and B is a real (BC) 2 CB anti- symmetric n x n matrix. Then 
the ground-state energy gap of TL{s) is bounded from below by a polynomial in 1/n and s. 

Proof. The proof is analogous to that of Corollary IIV. 31 □ 

V. REPRESENTATIONS OF THE HAMILTONIANS 

Using the Jordan- Wigncr transformation, we can define Hamiltonians using other kinds of particle operators that 
also have large ground-state energy gaps. In the Hubbard model of free electrons, interaction terms such as ejej. — CjC k 
do not occur because they do not conserve the number of electrons. However, they may occur in spin systems 
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transformed into fermionic representations. The best-known example is the Hamiltonian resulting from the Jordan- 
Wigner transformation applied to the XY-model Q. Let us first identify all the Hamiltonians that, under the 
Jordan- Wigncr transformation [26] 



= (-iy- 1 afal..aU 



v ' 



4 = (-iy-^...aU ] , (39) 



yield a quadratic fermionic Hamiltonian in the form of l[7|. Theorem IV. II is equivalent to the result in [271 . p. 4], but 
using a different basis representation. 

Theorem V.l (Quadratic fermionic Hamiltonians represented with Pauli operators). There is a bijection between 
Hamiltonians on n qubits of the form 

n 

U X W ^ + E ■ ■ ■ <%-i°i + E W k*°j°j+i ■ ■ ■ > ( 4 °) 

j=l k>j k>j 

where the coefficients Wj^k are real, and Hamiltonians of the form 

n 

n E A ^ { c h - c A) + B J,k (44 - , (41) 
i,k=i 

where {cj : j = 1, ...,n}, defined by i39\) . satisfy the FCRs, A is a real symmetric n x n matrix, and B is a real 
anti- symmetric matrix. The bijection is given by the invertible transformation 

A 3,3 = W j;j , 

(_l)m+l 

Aj,j+m = A h j +m = - (Wjj+ m + Wj +m j) , 

Bj,j+m = -Bj.j +m = [Wj,j+ m - Wj+mj) . (42) 

Proof. Apply (|3"9"|) to (j¥T|) . and use the commutation relations for Pauli operators to simplify the result. □ 

Using Theorem IV. 1| we can restate earlier results in a surprisingly simple manner. First, observe that application 
of Theorem IV. II to Ho defined in (JTSJ) yields 

n 

The ground state of 7io is the configuration with each particle in a spin-down eigenstate of <r z . If Hp is in the form 
of ([301), then so is the Hamiltonian evolution 

n 

H{s) = {1- s)Y^°-j +sH P , (44) 
j'=i 

for < s < 1. 

Next, observe that, up to sign, the elements of the matrix W are the same those of A + B. So to find the ground- 
state energy gap for a Hamiltonian that can be written in the form of P0"]) . we only need to apply the necessary sign 
changes to the elements of W, find the least non-zero singular value of the resulting matrix, and multiply by two. 
Thus Theorem IIII. 31 can be applied to the Hamiltonians in ([^0]) yielding a simple result: 

Theorem V.2 (Ground-state energy gaps of Hamiltonians defined using Pauli operators with Gaussian coefficients). 
Let H be defined by J^0[ ), where the elements of W are N(0,1) and independent. Let 7 be the ground-state energy gap 
ofH. Then, for large n, wy 2 /A converges in distribution to the probability density function 

p(x) = I±vJVW 2 +^> . (45) 
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Proof. Observe that the entries of W are, up to sign, those of A + B as defined by Theorem IV. 11 Thus A + B has 
independent N(0,1) entries, so we have the same proof as Theorem IIII. 31 □ 

The universal two-dimensional cluster state may be found in polynomial time using AQC (28j . The one-dimensional 
cluster state, while not universal for quantum computing, is useful for gaining intuition about cluster states p9j . The 
third-order interaction terms in (|4"U|) are exactly the stabilizers of the one-dimensional cluster state [3(| ■ In fact 

n-l 

H = - 2_, <7 3 CT 3+l (7 j+2 + {- 1 ) CT l f7 2 CT 3-- CT n-2<-l + (-l) ^2 a 3 a 4 ■ ■ ■ ( 46 ) 

is a Hamiltonian whose ground state is the one-dimensional cluster state. By Theorem IV. II A and B matrices 
corresponding to Ti. are circulant, and so by Theorem IIV.2I the one-dimensional cluster state can be realized in 
polynomial time. 

In general, we can define Fermi operators using spin-iS operators provided 2S +1 = 2™ for some n [l2| . Let us 
consider S = 3/2. Using n spin-3/2 particles, we can define 2n Fermi operators by 

(47) 

I) 2 ] , (48) 

where S x , S y , and S z are spin-3/2 operators and 5* = S x ± iS y . While the standard Jordan- Wigner transform 
applied to a one-dimensional chain of spin-1/2 particles results in a tridiagonal B matrix, the spin-3/2 transform 
applied to a one-dimensional chain of spin-3/2 particles yields a pcntadiagonal B matrix. 



f'2. 



1 

7s 



n 

2 

s? ) s: 



n 

k<j 



-0 



VI. CONCLUSION 



We showed that polynomially-large ground-state energy gaps are rare in many-body quantum Hermitian operators, 
but the gap is generically polynomially-large for quadratic fermionic Hamiltonians. We then proved that adiabatic 
quantum computing can realize the ground states of Hamiltonians with certain random interactions, as well as the 
ground states of one, two, and three-dimensional fermionic interaction lattices, in polynomial time. Finally, we used 
the Jordan- Wigner transformation to show that our results can be restated with Pauli operators in a surprisingly 
simple manner, and also consider a related spin-3/2 transformation. 

Since quadratic fermionic Hamiltonian evolutions are classically simulatable, the adiabatic quantum computations 
in Section UVl are simulatable. Thus we have provided a polynomial-time classical algorithm for finding properties of 
the ground states of certain random-interaction Hamiltonians and fermionic interaction lattices in one, two, and three 
dimensions. 

It should be noted that the Jordan- Wigner transformation can be generalized to higher dimensions, e.g. [3lj . 
Interesting results may follow from application of these alternate transformations to our theorems. 

Some fermionic systems may only approximately decouple into "non-interacting quasiparticles" , unlike the ex- 
act decouplings studied here. These systems may be "approximately" classically simulatable have "approximately" 
polynomially-large ground state energy gaps. This may be interesting to explore. 

Also, in principle, the decoupling-transformation may be classically difficult to find, whereas in the Hamiltonians 
studied here it is known how to find them. It is interesting whether similar results could be obtained for systems 
where the transformation is difficult to find explicitly. 



APPENDIX A: PROOF OF THE LIEB ET AL. THEOREM 



For completeness we include a proof of Theorem III. II We first need the property that the fermionic commutation 
relations are preserved under certain unitary transformations. 

Theorem A.l (Unitary transformations). Suppose the operators {cj : j = 1, . . . ,n} obey the FCRs. Let 

T=( U V V u). <-) 
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where U and V are real n x n matrices, and suppose T is unitary. Define the set of operators {r/j : j — 1, . . . , n} by 



where $A2\I denotes the transformation 



i 



\ r& J 



C2 



T 



V 4 / 



(A2) 



Tj,i c i 



T 



i=l 



Then {r/j : j = 1, . . . , n} also obey the FCRs. 

Proof. The proof follows from substituting the definitions of {r/j : j = 1 
commutation relations on {cj, ct : j = 1 . . .n}. 

Now we are ready to prove Theorem III. II 

Proof of Theorem I71.il We write Equation ([7]) as 



(A3) 



i} into the FCRs, and using the known 

□ 



C2 



H 



= (c\ ct 



ci c 2 



A B 

-B -A 



(A4) 



The theorem is equivalent to showing there are solutions to 



A B 

-B -A 



1({X + Y)(X-Y)\(K 0\ 1 
2\(X-Y) (X + Y) J I —A J 2 



(x + y) (x-y) 

(X-F) (A + F) 



(A5) 



for some non-negative real n x n diagonal matrix A, where X and Y are unitary. If so, then substituting Equation 
[5|) into Equation (|A4[) and using the definition of r)k, we get 



n 



(A6) 



fc=i 



Further, by Theorem lA.il {% : fc = 1 . . . n} satisfy the FCRs. So we can apply the FCRs to the second term in each 
summand to get Equation ((9|). 

Now we set about finding solutions to Equation (|A5j) . We rewrite it for convenience as: 



X + Y X-Y 
X -Y X + Y 



A B 

-B -A 



A 

-A 



X + Y X-Y 
X-Y X + Y 



(A7) 
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Equation (I A7|) is equivalent to the following four equations: 

(X + Y)A- (X -Y)B = A(X + Y) , (A8) 

(X + Y)B - (X - Y)A = A(X - Y) , (A9) 

(X — Y)A — (X + Y)B = —A(X — Y) , (A10) 

(X — Y)B — (X + Y)A = —A(X + Y) . (All) 

Evidently only two of the equations arc independent. Adding and subtracting Equations (|A8[) and (|A10[) yields 

X(A-B)=AY , (A12) 

Y(A + B)=AX. (A13) 

We can left-multiply by A to get 

AX (A — B) = A 2 Y , (A14) 

AY (A + B) = A 2 X , (A15) 

and then substitute Equation (|A13[) into Equation (|A14|) and Equation (|A12|) into Equation (|A15[) to get the pair of 
eigen-decomposition equations 

Y(A + B)(A-B) = A 2 Y , (A16) 
X(A — B){A + B) = A 2 X . (A17) 

Since A is real symmetric and B is real anti-symmetric, (A + B)^ = A — B and so (A — B)(A + B) and (A + B)(A — B) 
are symmetric positive semi-definite. So there is always a unitary X and Y with non-negative diagonal A 2 satisfying 
Equations (jA"l"6)) and (fAlT)) . □ 
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